ABSTRACT
Comprehensive analyses showed that SARS-CoV-2 infection caused COVID-19 and induced strong immune responses and sometimes severe illnesses. However, cellular features of recovered patients and long-term health consequences remain largely unexplored. In this study, we collected peripheral blood samples from nine recovered COVID-19 patients (median age of 36 years old) from Hubei province, China, 3 months after discharge as well as 5 age- and gender-matched healthy controls; and carried out RNA-seq and whole-genome bisulfite sequencing to identify hallmarks of recovered COVID-19 patients. Our analyses showed significant changes both in transcript abundance and DNA methylation of genes and transposable elements (TEs) in recovered COVID-19 patients. We identified 425 upregulated genes, 214 downregulated genes, and 18,516 differentially methylated regions (DMRs) in total. Aberrantly expressed genes and DMRs were found to be associated with immune responses and other related biological processes, implicating prolonged overreaction of the immune system in response to SARS-CoV-2 infection. Notably, a significant amount of TEs was aberrantly activated and their activation was positively correlated with COVID-19 severity. Moreover, differentially methylated TEs may regulate adjacent gene expression as regulatory elements. Those identified transcriptomic and epigenomic signatures define and drive the features of recovered COVID-19 patients, helping determine the risks of long COVID-19, and guiding clinical intervention.
ABSTRACT
BACKGROUND: The coronavirus disease 2019 (COVID-19) outbreak caused by severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) has become an ongoing pandemic. Understanding the respiratory immune microenvironment which is composed of multiple cell types, together with cell communication based on ligand-receptor interactions is important for developing vaccines, probing COVID-19 pathogenesis, and improving pandemic control measures. METHODS: A total of 102 consecutive hospitalized patients with confirmed COVID-19 were enrolled in this study. Clinical information, routine laboratory tests, and flow cytometry analysis data with different conditions were collected and assessed for predictive value in COVID-19 patients. Next, we analyzed public single-cell RNA-sequencing (scRNA-seq) data from bronchoalveolar lavage fluid, which offers the closest available view of immune cell heterogeneity as encountered in patients with varying severity of COVID-19. A weighting algorithm was used to calculate ligand-receptor interactions, revealing the communication potentially associated with outcomes across cell types. Finally, serum cytokines including IL6, IL1ß, IL10, CXCL10, TNFα, GALECTIN-1, and IGF1 derived from patients were measured. RESULTS: Of the 102 COVID-19 patients, 42 cases (41.2%) were categorized as severe. Multivariate logistic regression analysis demonstrated that AST, D-dimer, BUN, and WBC were considered as independent risk factors for the severity of COVID-19. T cell numbers including total T cells, CD4+ and CD8+ T cells in the severe disease group were significantly lower than those in the moderate disease group. The risk model containing the above mentioned inflammatory damage parameters, and the counts of T cells, with AUROCs ranged from 0.78 to 0.87. To investigate the molecular mechanism at the cellular level, we analyzed the published scRNA-seq data and found that macrophages displayed specific functional diversity after SARS-Cov-2 infection, and the metabolic pathway activities in the identified macrophage subtypes were influenced by hypoxia status. Importantly, we described ligand-receptor interactions that are related to COVID-19 serverity involving macrophages and T cell subsets by communication analysis. CONCLUSIONS: Our study showed that macrophages driving ligand-receptor crosstalk contributed to the reduction and exhaustion of CD8+ T cells. The identified crucial cytokine panel, including IL6, IL1ß, IL10, CXCL10, IGF1, and GALECTIN-1, may offer the selective targets to improve the efficacy of COVID-19 therapy. TRIAL REGISTRATION: This is a retrospective observational study without a trial registration number. Video Abstract.